Estudo fenomenológico e numérico do escoamento estratificado óleo-água ondulado e com mistura na interface

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura.

Modelagem do absorvedor e do gerador de ciclos de refrigeração por absorção de calor com o par amônia/água baseados na tecnologia de filme descendente sobre placas inclinadas.

Esse trabalho constitui o desenvolvimento da modelagem térmica e simulação por métodos numéricos de dois componentes fundamentais do ciclo de refrigeração por absorção de calor com o par amônia/água: o absorvedor e o gerador. A função do absorvedor é produzir mistura líquida com alta fração mássica de amônia a partir de mistura líquida com baixa fração mássica de amônia e mistura vapor mediante retirada de calor. A função do gerador é produzir mistura líquido/vapor a partir de mistura líquida mediante o fornecimento de calor. É proposto o uso da tecnologia de filmes descendentes sobre placas inclinadas e o método de diferenças finitas para dividir o comprimento da placa em volumes de controle discretos e realizar os balanços de massa, espécie de amônia e energia juntamente com as equações de transferência de calor e massa para o filme descendente. O objetivo desse trabalho é obter um modelo matemático simplificado para ser utilizado em controle e otimização. Esse modelo foi utilizado para calcular as trocas de calor e massa no absorvedor e gerador para diversas condições a partir de dados operacionais, tais como: dimensões desses componentes, ângulo de inclinação da placa, temperatura de superfície e condições de entrada da fase líquida e vapor. Esses resultados foram utilizados para estabelecer relações de causa e efeito entre as variáveis e parâmetros do problema. Os resultados mostraram que o ângulo de inclinação da placa ótimo tanto para o absorvedor como para o gerador é a posição vertical, ou 90°. A posição vertical proporciona o menor comprimento de equilíbrio (0,85 m para o absorvedor e 1,27 m para o gerador com as condições testadas) e se mostrou estável, pois até 75° não foram verificadas variações no funcionamento do absorvedor e gerador. Dentre as condições testadas para uma placa de 0,5 m verificou-se que as maiores efetividades térmicas no absorvedor e gerador foram respectivamente 0,9 e 0,7 e as maiores efetividades mássicas no absorvedor e gerador foram respectivamente 0,6 e 0,5. É esperado que os dados obtidos sejam utilizados em trabalhos futuros para a construção de um protótipo laboratorial e na validação do modelo.

Estudo de transistores de tunelamento induzido por efeito de campo (TFET) construídos em nanofio.

Esse trabalho de mestrado teve como estudo o transistor Túnel-FET (TFET) fabricado em estrutura de nanofio de silício. Este estudo foi feito de forma teórica (simulação numérica) e experimental. Foram estudadas as principais características digitais e analógicas do dispositivo e seu potencial para uso em circuitos integrados avançados para a próxima década. A análise foi feita através da extração experimental e estudo dos principais parâmetros do dispositivo, tais como inclinação de sublimiar, transcondutância (gm), condutância de saída (gd), ganho intrínseco de tensão (AV) e eficiência do transistor. As medidas experimentais foram comparadas com os resultados obtidos pela simulação. Através do uso de diferentes parâmetros de ajuste e modelos de simulação, justificou-se o comportamento do dispositivo observado experimentalmente. Durante a execução deste trabalho estudou-se a influência da escolha do material de fonte no desempenho do dispositivo, bem como o impacto do diâmetro do nanofio nos principais parâmetros analógicos do transistor. Os dispositivos compostos por fonte de SiGe apresentaram valores maiores de gm e gd do que aqueles compostos por fonte de silício. A diferença percentual entre os valores de transcondutância para os diferentes materiais de fonte variou de 43% a 96%, sendo dependente do método utilizado para comparação, e a diferença percentual entre os valores de condutância de saída variou de 38% a 91%. Observou-se também uma degradação no valor de AV com a redução do diâmetro do nanofio. O ganho calculado a partir das medidas experimentais para o dispositivo com diâmetro de 50 nm é aproximadamente 45% menor do que o correspondente ao diâmetro de 110 nm. Adicionalmente estudou-se o impacto do diâmetro considerando diferentes polarizações de porta (VG) e concluiu-se que os TFETs apresentam melhor desempenho para baixos valores de VG (houve uma redução de aproximadamente 88% no valor de AV com o aumento da tensão de porta de 1,25 V para 1,9 V). A sobreposição entre porta e fonte e o perfil de dopantes na junção de tunelamento também foram analisados a fim de compreender qual combinação dessas características resultariam em um melhor desempenho do dispositivo. Observou-se que os melhores resultados estavam associados a um alinhamento entre o eletrodo de porta e a junção entre fonte e canal e a um perfil abrupto de dopantes na junção. Por fim comparou-se a tecnologia MOS com o TFET, obtendo-se como resultado um maior valor de AV (maior do que 40 dB) para o TFET.

Estudo do processo de formação do cavaco durante o torneamento e sua relação com a microestrutura utilizando o método dos elementos finitos.

O objetivo desse trabalho foi avaliar o processo de formação do cavaco durante o torneamento utilizando simulação numérica pelo método dos elementos finitos. Para realizar o estudo foram definidos dois tipos de aços inoxidáveis austeníticos, um com matriz metálica sem a presença significativa de inclusões, do tipo ABNT 304, e outro com a presença de inclusões não metálicas, do tipo ABNT 303. O estudo foi focado nos mecanismos de formação e ruptura do cavaco, na determinação das forças de usinagem, no campo de tensões, deformações, e temperaturas durante o processo, que foram relacionados com aspectos e características da microestrutura do material. Os resultados obtidos foram comparados com as forças de usinagem experimentais, com a espessura e morfologia do cavaco. O desenvolvimento do trabalho, de acordo com a metodologia adotada, foi realizado em diferentes etapas. Inicialmente foi elaborado e aplicado um modelo de simulação da usinagem considerando o material homogêneo. Em outra etapa, foi realizada a modelagem de uma microestrutura submetida a um estado de tensão e deformação semelhante ao encontrado na simulação da usinagem realizada com material homogêneo. Os resultados mostraram que as partículas das inclusões maiores, alongadas, e em maior quantidade aumentam a tensão e a deformação na microestrutura. As elevadas temperaturas obtidas na usinagem dos aços inoxidáveis austeníticos aumentam a ductilidade dos sulfetos, esses se deformam em compressão junto com a matriz, e têm um efeito limitado como agente de redução dos esforços de usinagem. Por outro lado, os sulfetos facilitam a etapa de ruptura do cavaco em tensões trativas, e tendem a se romper facilitando o processo de quebra.

Models and methods for the direct simulation of rough and micropatterned surfaces

Friction in hydrodynamic bearings are a major source of losses in car engines ([69]). The extreme loading conditions in those bearings lead to contact between the matching surfaces. In such conditions not only the overall geometry of the bearing is relevant, but also the small-scale topography of the surface determines the bearing performance. The possibility of shaping the surface of lubricated bearings down to the micrometer ([57]) opened the question of whether friction can be reduced by mean of micro-textures, with mixed results. This work focuses in the development of efficient numerical methods to solve thin film (lubrication) problems down to the roughness scale of measured surfaces. Due to the high velocities and the convergent-divergent geometries of hydrodynamic bearings, cavitation takes place. To treat cavitation in the lubrication problem the Elrod- Adams model is used, a mass-conserving model which has proven in careful numerical ([12]) and experimental ([119]) tests to be essential to obtain physically meaningful results. Another relevant aspect of the modeling is that the bearing inertial effects are considered, which is necessary to correctly simulate moving textures. As an application, the effects of micro-texturing the moving surface of the bearing were studied. Realistic values are assumed for the physical parameters defining the problems. Extensive fundamental studies were carried out in the hydrodynamic lubrication regime. Mesh-converged simulations considering the topography of real measured surfaces were also run, and the validity of the lubrication approximation was assessed for such rough surfaces.

Comportamento tribo-mecânico e desgaste adesivo de materiais em nanoescala: análises por dinâmica molecular e mecânica do contínuo.

As formulações baseadas na mecânica do contínuo, embora precisas até certo ponto, por vezes não podem ser utilizadas, ou não são conceitualmente corretas para o entendimento de fenômenos em escalas reduzidas. Estas limitações podem aparecer no estudo dos fenômenos tribológicos em escala nanométrica, que passam a necessitar de novos métodos experimentais, teóricos e computacionais que permitam explorar estes fenômenos com a resolução necessária. Simulações atomísticas são capazes de descrever fenômenos em pequena escala, porém, o número necessário de átomos modelados e, portanto, o custo computacional - geralmente torna-se bastante elevado. Por outro lado, os métodos de simulação associados à mecânica do contínuo são mais interessantes em relação ao custo computacional, mas não são precisos na escala atômica. A combinação entre essas duas abordagens pode, então, permitir uma compreensão mais realista dos fenômenos da tribologia. Neste trabalho, discutem-se os conceitos básicos e modelos de atrito em escala atômica e apresentam-se estudos, por meio de simulação numérica, para a análise e compreensão dos mecanismos de atrito e desgaste no contato entre materiais. O problema é abordado em diferentes escalas, e propõe-se uma abordagem conjunta entre a Mecânica do Contínuo e a Dinâmica Molecular. Para tanto, foram executadas simulações numéricas, com complexidade crescente, do contato entre superfícies, partindo-se de um primeiro modelo que simula o efeito de defeitos cristalinos no fenômeno de escorregamento puro, considerando a Dinâmica Molecular. Posteriormente, inseriu-se, nos modelos da mecânica do contínuo, considerações sobre o fenômeno de adesão. A validação dos resultados é feita pela comparação entre as duas abordagens e com a literatura.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Modelagem e análise de desempenho de sistema para geração de energia elétrica através de ondas marítimas.

Mediante a crescente necessidade de aumento na oferta de energia elétrica devido à constante elevação na demanda mundial, esta dissertação avalia o desempenho de um sistema conversor de energia de ondas marítimas em energia elétrica. O sistema em análise é o de coluna de água oscilante com turbina de dupla ação instalado na costa. Utiliza-se um modelo regular de ondas como perturbação à dinâmica de uma câmara semi-submersa gerando fluxo de ar através de uma turbina à ar de dupla ação. O sistema final é não linear e com parâmetros variantes no tempo. A dissertação investiga possibilidades para o aumento do rendimento da turbina em diferentes condições de mar através do método de simulação numérica. Após a modelagem física e matemática do sistema escolhido, inicia-se a síntese de um controlador proporcional derivativo para controle da pressão de ar na turbina em torno da pressão ideal de trabalho da mesma. A análise inclui o comparativo entre os resultados do sistema com e sem controlador e a avaliação de robustez utilizando ondas com amplitude variável. O trabalho apresenta ainda propostas de otimização do sistema para trabalhar em condições similares a região de Pecém no Brasil. Pelos resultados obtidos nas simulações, conclui-se que o rendimento e a robustez do sistema podem melhorar utilizando um sistema controlado. O rendimento do sistema poderá ainda ser otimizado para a região de instalação.

Ensino-aprendizagem de enfermagem em simulação clínica: desenvolvendo competência profissional para prevenção de úlceras por pressão

O desenvolvimento de competências profissionais em processos de ensino-aprendizagem insere a simulação na perspectiva tecnológica realística de situações clínicas para atingir objetivos educacionais, visando segurança e qualidade. A ação competente operacionaliza saberes pela combinação de conhecimentos, habilidades e atitudes. A competência avaliação de risco para úlcera por pressão é privativa do enfermeiro, o agravo é indicador de qualidade e segurança do paciente. O objetivo do estudo foi avaliar o processo de ensino-aprendizagem ante a estratégia de simulação clínica, visando o desenvolvimento da competência profissional avaliação de risco para úlceras por pressão e seus respectivos conhecimentos, habilidades e atitudes. Fundamentou-se, teoricamente, o estudo, nos preceitos da andragogia, no modelo de Simulação da National League for Nursing/Jeffries, nas concepções de competências de Le Boterf e recomendações para prevenção de úlcera por pressão do NPUAP/ EPUAP. Trata-se de estudo descritivo com análise quanti e qualitativa, desenvolvido no Laboratório de Simulação da Universidade Federal do Piauí, em Teresina, com estudantes de enfermagem do último ano do curso. Constou de três etapas: construção e validação de conteúdo do instrumento para conhecimentos, habilidades e atitudes da competência, elaboração do plano de aula e sequência didática do cenário de simulação, aplicação de pré/pós-teste sobre o componente \"conhecimento\", aula expositiva, execução do cenário de simulação pelos estudantes, avaliação dos cenários pelo comitê de juízes e debriefing por grupo focal, após o cenário. Verificou-se confiabilidade do construto e do instrumento de medida da competência com concordância predominando de 80 a 100% de pertinência. O instrumento constou de 32 itens: conhecimento (14), habilidades (8) e atitudes (10), com escala de resposta de cinco pontos, de nada (1) a extremamente (5). Inicialmente o estudo contou com 35 estudantes. Desses, 29 (82,8%) participaram da simulação. Predominou o sexo feminino 22 (76%), faixa etária de 21 a 24 anos, 24 (82,7%) e procedentes de Universidades Públicas do Estado do Piauí (83%). Os resultados mostraram melhores níveis de combinação de saberes para o componente conhecimentos, após a aplicação das estratégias de ensino autorreferidos pelos estudantes. Na avaliação dos conhecimentos, habilidade e atitudes, durante o cenário, houve discordância entre os juízes: para o juiz 1 predominaram parâmetros negativos da escala (nada e muito pouco) e para o juiz 2 e 3 os positivos (bastante e extremamente). As comparações pelo teste de variância ANOVA mostraram diferenças significativas entre juiz 1 e 2 (p-valor = < 0,01) e juiz 1 e 3 (p-valor=<0,01) para todos os componentes. A análise dos dados obtidos, durante o debriefing, originou cinco categorias com suas respectivas unidades de registro temáticas pela análise de conteúdo de Bardin. Os resultados sugerem que a estratégia é capaz de resgatar o raciocínio operativo dos estudantes durante a ação, desenvolver pensamento crítico-reflexivo sobre a competência, identificar lacunas de aprendizagem, promover satisfação aos estudantes e melhorar a autoimagem profissional. Conclui-se que a estratégia possibilita o desenvolvimento da competência avaliação de risco para úlcera por pressão nas dimensões de saberes (conhecimento), fazeres (habilidades) e querer-agir, saber-agir e poder-agir (atitudes), explorados neste estudo.

Julgamento clínico e raciocínio diagnóstico de estudantes de enfermagem em simulação clínica de alta-fidelidade

O ensino dos processos de julgamento clínico e de raciocínio diagnóstico para estudantes de enfermagem torna-se cada vez mais importante para a qualificação dos cursos de graduação em enfermagem e alcance do compromisso com a formação clínica por excelência do enfermeiro, em uma realidade de saúde cada vez mais complexa. Os objetivos deste estudo foram identificar o julgamento clínico e o raciocínio diagnóstico de estudantes de enfermagem, correlacionar estes dois processos e identificar diferenças entre estudantes da fase intermediária e os concluintes do curso, de uma Escola do interior paulista. Para avaliar o julgamento clínico, construímos um cenário de simulação clínica de alta-fidelidade, representando uma paciente com anemia falciforme em crise de dor e, ainda, traduzimos e adaptamos à cultura brasileira o instrumento Lasater Clinical Judgment Rubric (LASATER, 2007); para avaliar o raciocínio diagnóstico, utilizamos o instrumento já adaptado por Rodrigues (2012) denominado de Inventário de Raciocínio Diagnóstico (BORDAGE; GRANT; MARSDEN, 1990). Os resultados demonstraram que os estudantes de enfermagem apresentaram, predominantemente, nível Proficiente na maioria das dimensões de julgamento clínico (66,7% dos estudantes do grupo concluinte e 56,5% dos estudantes do grupo intermediário). Já para o raciocínio diagnóstico, a maioria dos estudantes foi considerada com ampla habilidade para realizar diagnósticos de enfermagem (91,3% dos estudantes do grupo intermediário e 83,4% dos estudantes do grupo concluinte). Destaca-se que dos concluintes 11,1% apresentaram habilidade máxima. Os estudantes do último ano de graduação em enfermagem apresentaram desempenho superior na fase de interpretação do julgamento clínico (p=0,021). Não se observou diferença entre os grupos para o raciocínio diagnóstico (p=0,334). Houve moderada correlação entre julgamento clínico e raciocínio diagnóstico; e ainda, a fase de reconhecimento do julgamento clínico apresentou-se moderadamente correlacionada ao processo de raciocínio diagnóstico. Considerando que o raciocínio diagnóstico está presente no processo de julgamento clínico, principalmente no momento da investigação do caso clínico (fase Reconhecimento do julgamento clínico) e que as habilidades de raciocinio diagnóstico se manifestam predominantemente nesta fase, a compreensão e o desenvolvimento destes processos pelos estudantes devem ser valorizados nos currículos de graduação em enfermagem

Estudo dos parâmetros dosimétricos de sementes de Iodo-125 desenvolvidas pelo IPEN-CNEN/SP utilizadas em braquiterapia por simulação computacional pelo método de Monte Carlo

As expectativas da Organização Mundial de Saúde para o ano de 2030 são que o número de mortes por câncer seja de aproximadamente 13,2 milhões, evidenciando a elevada parcela desta doença no problema de saúde mundial. Com relação ao câncer de próstata, de acordo com o Instituto Nacional do Câncer, o número de casos diagnosticados no mundo em 2012 foi de aproximadamente 1,1 milhão, enquanto que no Brasil os dados indicam a incidência de 68 mil novos casos. O tratamento deste tipo de neoplasia pode ser realizado com cirurgia (prostatectomia) ou radioterapia. Dentre a radioterapia, podemos destacar a técnica de braquiterapia, a qual consiste na introdução (implante) de pequenas fontes radioativas (sementes) no interior da próstata, onde será entregue um valor elevado de dose no volume de tratamento e baixa dose nos tecidos ao redor. No Brasil, a classe médica estima uma demanda de aproximadamente 8000 sementes/mês, sendo o custo unitário de cada semente de pelo menos U$ 26,00. A Associação Americana de Físicos na Medicina publicou alguns documentos descrevendo quais parâmetros e análises devem ser realizadas para avaliações da distribuição de dose, como por exemplo, os parâmetros Constante de taxa de dose, Função radial e Função de anisotropia. Estes parâmetros podem ser obtidos através de medidas experimentais da distribuição de dose ou por simulações computacionais. Neste trabalho foram determinados os parâmetros dosimétricos da semente OncoSeed-6711 da empresa Oncura-GEHealthcare e da semente desenvolvida pelo Grupo de Dosimetria de Fontes de Braquiterapia do Centro de Tecnologia das Radiações (CTR IPEN-CNEN/SP) por simulação computacional da distribuição de dose utilizando o código MCNP5, baseado no Método de Monte Carlo. A semente 6711 foi modelada, assim como um sistema dosimétrico constituído por um objeto simulador cúbico de 30x30x30 cm3 preenchido com água. Os valores obtidos da semente 6711 foram comparados com alguns apresentados na literatura, onde o parâmetro Constante de taxa de dose apresentou erro relativo em relação ao valor publicado no TG- 43 de 0,1%, sendo que os outros parâmetros analisados também apresentaram boa concordância com os valores publicados na literatura. Deste modo, pode-se considerar que os parâmetros utilizados nas simulações (espectro, modelagem geométrica e avaliação de resultados) estão compatíveis com outros estudos, sendo estes parâmetros também utilizados nas simulações da semente do IPEN. Considerando as análises de incerteza estatística, os valores obtidos da semente do IPEN são semelhantes aos valores da semente 6711.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Avaliação de conflitos de usos múltiplos de água na bacia hidrográfica do Rio Jaguari - SP, através de simulação computacional

No gerenciamento de recursos hídricos, a bacia hidrográfica é considerada a unidade natural de gestão. Suas características físicas definem um espaço onde a água converge para uma saída comum. O desenvolvimento integrado de uma bacia, bem como os possíveis conflitos ocasionados por demandas concorrentes do mesmo recurso, devem ser administrados respeitando esse conceito de unidade. No caso específico da bacia hidrográfica do Rio Jaguari, existe uma crescente preocupação relacionada aos usos múltiplos da água superficial. O rio Jaguari, importante afluente do Paraíba do Sul no estado de São Paulo, tem suas águas represadas para geração de energia elétrica e regularização de vazões. Nos últimos anos, tem-se constatado um crescente rebaixamento dos níveis operacionais da represa do Jaguari. Segundo estudos da sociedade civil organizada local, a tendência é de esgotamento do reservatório em poucos anos. Este trabalho aborda a questão do rebaixamento dos níveis de Jaguari através da aplicação do software DHI Mike Basin 2000. Trata-se de um simulador genérico para sistemas de recursos hídricos, de ampla aplicação. A bacia do Jaguari é representada em um modelo matemático e são simulados quatro cenários distintos de usos de água na bacia. Verifica-se que o problema é real e requer uma intervenção multi-institucional para ser solucionado.

Análise de conflitos de usos múltiplos da água auxiliada por simulação computacional: bacia do rio Sapucaí - Mirim/Grande

Na mesma proporção que fontes de água potável são reduzidas, a competição por elas tem aumentado, conseqüência imediata da expansão de áreas urbanas e de novas demandas agro-industriais. Os conflitos gerados por esta competição tendem a comprometer o desenvolvimento sustentável e o progresso de uma forma geral, uma vez que resultam em perda de eficácia para vultuosos investimentos públicos e privados e em prejuízos para os usuários menos estruturados. Considerando os fatos acima, o objeto deste trabalho é o estabelecimento de cenários de alocação das disponibilidades hídricas entre usuários de água da bacia do Sapucaí-Mirim/Grande, visando determinar a convivência harmoniosa das várias finalidades de uso provenientes do mesmo recurso hídrico. Tal meta poderá ser alcançada através da análise dos indicadores de desempenho - confiabilidade, vulnerabilidade e resiliência - do sistema os quais serão obtidos com o auxílio do modelo de simulação IRAS (Interactive River Aquifer Simulation).

Estudo da solidificação equiaxial utilizando o modelo do campo de fases tridimensional.

Este trabalho apresenta um estudo da solidificação de metais puros utilizando o modelo de campo de fases. O modelo é utilizado para simular a solidificação com o intuito de obter a morfologia da interface sólido-líquido sob diversas condições de transferência de calor. Foram realizados testes de validação comparando as morfologias da interface sólido-líquido obtida com as morfologias apresentadas em trabalhos anteriores para os casos bi e tridimensionais. O modelo do campo de fases adotado consiste principalmente de duas equações diferenciais: uma para calcular a variável de campo de fases e outra para calcular o campo de temperaturas. As equações foram solucionadas numericamente para um oitavo do domínio devido a simetria do problema. Os cálculos do modelo indicam que um sólido esférico com um raio inicial menor que o raio crítico de nucleação refunde. Entretanto uma esfera de raio maior cresce. Quando o sólido inicial cresce em uma malha numérica relativamente grosseira, a forma do sólido desvia da forma esférica devido perturbações na interface sólido-líquido. Quando a malha é refinada, as perturbações não são detectadas; contudo, quando introduzidas artificialmente as perturbações crescem e distorcem o formato esférico.